Constructions of Quantum Convolutional Codes 



Markus Grassl 
Institut fur Algorithmen und Kognitive Systeme 
Fakultat fur Informatik, Universitat Karlsruhe (TH) 
Am Fasanengarten 5, 76 128 Karlsruhe, Germany 
Email: grassl@ira.uka.de 



Martin Rotteler 
NEC Laboratories America, Inc. 
4 Independence Way, Suite 200 

Princeton, NJ 08540, USA 
Email: mroetteler@nec-labs.com 



o 
o 

(N 



> : 

(N . 
00 ■ 

7— I 1 

m : 
o . 
r- ■ 
o ■ 

^ : 

i 

^— > ■ 

G ■ 
cd 

3 : 



Abstract — We address the problems of constructing quantum 
convolutional codes (QCCs) and of encoding them. The first 
construction is a CSS-type construction which allows us to find 
QCCs of rate 2/4. The second construction yields a quantum 
convolutional code by applying a product code construction 
to an arbitrary classical convolutional code and an arbitrary 
quantum block code. We show that the resulting codes have 
highly structured and efficient encoders. Furthermore, we show 
that the resulting quantum circuits have finite depth, independent 
of the lengths of the input stream, and show that this depth is 
polynomial in the degree and frame size of the code. 

I. Introduction 

Similar to the classical case a quantum convolutional code 
(QCC) encodes an incoming stream of quantum information 
into an outgoing stream. A basic theory of quantum convolu- 
tional codes obtained from infinite stabilizer matrices has been 
developed recently, see [13]. 

Only few constructions of quantum convolutional codes are 
known, see [2], [3], [4], [5], [6], [8], [13]. In this paper, we 
construct some new quantum convolutional codes using a CSS- 
type construction which uses the same principle as the CSS 
construction for block codes [12]. Furthermore, we revisit the 
product code construction introduced in [8] and show that for 
these codes the algorithm presented in [9] for computing a 
non-catastrophic encoder takes a particularly simple form. This 
allows us to show that the depth of the encoding circuit is 
polynomial in the frame size and the constraint length of the 
code. 

II. Quantum convolutional codes 

A. Basic definitions 

QCCs are defined as infinite versions of quantum stabilizer 
codes. The appropriate generalization of stabilizer block codes 
to QCCs is provided by the polynomial formalism introduced 
in [13]. We briefly sketch this approach^ 

The code is specified by its stabilizer which is a subgroup 
of the infinite version Qoo of the Pauli group, which consists 
of tensor products of generalized Pauli matrices acting on an 
semi-infinite stream of qudits. The stabilizer can be described 
by a matrix with polynomial entries 

5(D) = (X(D)\Z(D)) e ¥ q [D] {n - k ^ 2n . (1) 

'We describe the approach for q dimensional subsystems (qudits) which is 
a straightforward generalization of the binary case. 



Definition 1: Let C be a QCC defined by a full -rank sta- 
bilizer matrix as in eq. (Q~|). Then n is called the frame 
size, k the number of logical qudits per frame, and k/n 
the rate of the QCC. The constraint lengths are defined as 
vi = maxi<j<„(max(degXy(D),degZy (D))), the overall 
constraint length is the defined as the sum v — Y17=i v i> an d 
the memory m is given by m = maxi<;<„_fc vi. 

Like in the classical case, a QCC can also be described in 
terms of a semi-infinite stabilizer matrix S which has entries 
in F q x W q . First, we write 5(D) = J2?=o G i Dl usin g m + 1 
matrices Gq,G\, . . . , G m which are of size (n — k) x n each. 
Then we define the semi-infinite matrix 



5 := 
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(2) 



Note that 5 has a block band structure where each block is of 
size (n — k) x (m + l)n. A useful property of 5 is that every 
qudit in the semi-infinite stream of qudits is acted upon non- 
trivially by only a finite number of generators. Moreover, those 
generators have bounded support. Hence their eigenvalues 
can be measured as soon as the corresponding qudits have 
been received. Therefore, it is possible to compute the error 
syndrome for the quantum convolutional code online. 

There is a condition to check whether 5 is well-defined, i. e., 
if it defines a commutative subgroup of [13]. If S(D) = 
(X(D)\Z(D)) as in eq. ([TJ, then the condition of symplectic 
orthogonality of S translates to 



X{D)Z{1/Df - Z{D)X(l/D) t = 0. 



(3) 



Example 2: As an example we consider the QCC defined 
by the stabilizer matrix (see [5]) 



5(D) = 



1 + D 1 
D 



l + D 
D 





l + D 



D D 
l + D 1 



This code is derived from the classical Fzi-linear code gener- 
ated by (1 + D, 1 + cjD, 1 + uj 2 D). We can easily check self- 
orthogonality by computing X{D)Z{1/D) 1 - Z{D)X(1/D) t 
which turns out to be the 2 x 2 all zero matrix. Hence 
the code indeed is self-orthogonal to all shifted versions of 
itself, i. e., it defines a QCC where n = 3, k = 1, and 
in = 1. To illustrate the structure in terms of Pauli matrices 



TABLE I 

Action of various generalized Clifford operations. 
Conjugation by the unitary U corresponds to the 

ACTION OF U ON THE COLUMNS OF S(D) = (X(D)\Z(D)). 



unitary gate U 
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Fourier transform DFT 
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we consider the corresponding semi-infinite stabilizer matrix, 
which is given (in Pauli form) as follows: 



( X X X X Z Y 

Z Z Z Z Y X 

X X X X Z Y 

Z Z Z Z Y X 



\ 



X X X X Z Y 
Z Z Z Z Y X 



V 



/ 



It is easy to see that the QCC corresponding to S can correct 
an arbitrary number of errors, as long as they do not occur 
in bursts, meaning in this example that at least six unaffected 
qubits are between two erroneous ones. 

B. Encoding circuits 

In [7] it has been shown that for any block quantum error- 
correcting code C = [n, k, d\ q there is quantum circuit of 
polynomial size for encoding. In order to encode k qudits into 
n qudits, the circuit acts on the k input qudits and n— k ancillae 
which are initialized in the state |0). The input state can be 
described by a Z-only stabilizer matrix Sq — (X\Z) = (0|/0), 
where / is an (n—k) x (n—k) identity matrix. The operation of 
the encoding circuit corresponds to a transformation changing 
So into the stabilizer S of the quantum code. This idea can 
be adapted to quantum convolutional codes (see [9] for qubit 
codes). The encoding circuit can be realized by generalized 
Clifford gates whose action is summarized in Table H] for the 
gates and the corresponding actions, see [7, Theorem 2]. 



III. Efficient encoders for CSS type QCCs 

The CSS-like construction of QCCs uses two classical 
convolutional codes C\ = (n, k\) and C 2 = (n,n — k 2 ) with 
equal frame length n and C 2 C C\. The stabilizer matrix 
(X(D)\Z(D)) is of block diagonal form, given by 

H 2 (B) 
H 1 (D) 



e F q [D 



(n-fe+b)x2n 



(4) 



where Hi(D), H 2 (B) denote parity check matrices of C\ and 
C 2 , respectively. We assume that both Hi(D) and H 2 (B) 
correspond to non-catastrophic, delay-free encoders and have 
full ranks n — ki and k2, respectively. This implies that their 
Smith normal form is (1 0) with a suitable unit matrix / (see 
[10, Chapter 2]). In particular there are unimodular matrices 
Ai(D) G ¥ q [D] k2Xk2 and Bi(D) e 



?nxn such that 



A 1 (D)H 2 (D)B 1 (D) = (10). 



(5) 



There is an algorithm for computing the Smith normal form 
and the transformation matrices Ai(D) and Bi(D) whose bit- 
complexity is polynomial in the size and degree of the matrix 
B-z(D) [11]. This implies that the corresponding quantum 
circuit implementing the matrix B\(D) can be realized using 
polynomially many generalized Clifford gates. The trans- 
formed stabilizer matrix is given by 



(X'(D)\Z'(D)) = 





Z 2 (D) 



(6) 



The action of the gates results in a modified Z-part as well. 
From condition (0 it follows that Z^iB) is of the form 
Z 2 (B) = (0 Z' 2 {B)) with Z' 2 (D) e F j£)](«-fci)x(™-fc2). 
Using Fourier transformation (DFT) gates on the last n — k 2 
qudits, the stabilizer matrix reads 



(X'{D)\Z\D)) = 





Z' 2 {B) 



(7) 



Computing the Smith form of Z 2 (B) yields matrices A 2 (B) 
and B 2 {D) with A 2 (D)Z' 2 {D)B 2 {D) = (10). Another 
Fourier transform on the first n — k\ + k 2 qudits yields an 
Z-only stabilizer matrix of the form (0|/0). The resulting 
quantum code encodes ki — k 2 qudits per frame of size n. 
Overall, we obtain the following result. 

Theorem 3: Let C be a quantum convolutional code con- 
structed using the CSS-like construction from two classical 
convolutional codes C\ and C 2 with stabilizer matrix as in 
eq. (0). Denote the frame size with n and the constraint 
length with v. Then C has an encoding circuit whose depth is 
finite, i, e., does not depend on the length of the input stream. 
Furthermore, the depth of this circuit is upper bounded by 
poly(n, v). 

IV. CSS-TYPE QCCS OF RATE 2/4 

In [6], optimal quantum convolutional codes of rate 1/3 are 
listed which are based on self-orthogonal binary convolutional 
codes of rate 1/3. In order to construct quantum convolutional 
codes of rate 2 /4, we search for self-orthogonal binary convo- 
lutional codes C of rate 1/4 which have a dual code C7 X with 



2 



V 


9i(D) 


92(D) 


93(D) 


94(D) 


d 1 - 


N d ± 


3 


1100 


1110 


1001 


1101 


3 


2 


4 


11001 


11101 


10011 


10111 


4 


1 


4 


10001 


10101 


11011 


11111 


4 


1 


5 


110010 


111010 


100001 


110111 


5 


14 


6 


1010010 


1111010 


1000101 


1100111 


6 


63 


7 


10101001 


11111001 


10000011 


11000111 


6 


8 


8 


101100001 


100100101 


111110011 


111011011 


6 


2 


9 


1001001001 


1100111101 


1110110111 


1010101111 


7 


10 


10 


11011011001 


10100001101 


10011000011 


11001001111 


8 


67 


11 


101101100011 


101010110011 


111101001011 


110001101111 


8 


25 



Fig. 1 . Generators for self-orthogonal binary convolutional codes of rate 1/4 yielding quantum convolutional codes of rate 2/4 found by random search. 



high minimum distance d^. Applying the CSS construction 
with C\ — C 2 = C , we then obtain a quantum convolutional 
code of rate 2/4 and minimum distance d^~. 

The results of a randomized search for such codes is 
presented in Table 03 The entries of the generator matrix 
g(D) = (gi(D),g 2 (D),g 3 (D),g 4 (D)) of the code C are 
given in abbreviated form, listing the coefficients in increasing 
order. For example, the generator matrix of the first code with 
constraint length v = 3 is g(D) = (1 + D, 1 + D + D 2 , 1 + 
D 3 , 1 + D + D 3 ). The last column lists the number N d ± of 
sequences of minimum weight. Note that it is desirable to 
have as few sequences of minimum weight as possible. The 
size of the search space grows with 0(2 4l/ ), so we have only 
performed an exhaustive search up to constraint length v = 6, 
and a randomized search for larger values of v. 

V. Efficient encoders for product codes 

A. Product code construction 

The following theorem, taken from [8], allows to construct 
a quantum convolutional code using a classical convolutional 
code and a quantum code. 

Theorem 4: Let C\ = (raj., ki) p be a classical convolutional 
code over ¥ p with dual distance d^ and let G\(D) be a 
generator matrix of C\ corresponding to a non-catastrophic, 
delay-free encoder. Furthermore, let C be a quantum error- 
correcting code for q-dimensional quantum systems (q — pr) 
with minimum distance d 2 and stabilizer matrix S 2 = (X\Z) 
if C is a block code or S 2 = (X(D)\Z(D)) if C is a 
convolutional code. Then the stabilizer matrix 



G(D)=G 1 (D) ® P S 2 



(8) 



defines a quantum convolutional code with minimum distance 
d < min((ij L , d 2 ). 

The tensor product ® p corresponds to the Kronecker prod- 
uct of the stabilizer matrices. We use the index p to stress that 
the coefficients of the polynomials in the matrix G\(D) are 
in the prime field F p while the stabilizer matrix S 2 might be 
defined over an extension field ¥ q = F p * . 



B. Encoding product codes 

Instead of applying the general algorithm of [9] to the matrix 
G(D) in order to compute an encoding circuit for the product 
code, we will exploit the additional structure of the stabilizer 
matrix. The first step is to compute an inverse encoding circuit 
for the quantum code C with stabilizer S2. The quantum circuit 
corresponds to a symplectic transformation yielding the trivial 
Z-only stabilizer Sq = (0|/0). Note that the trivial stabilizer 
is of this form, regardless whether the code C is a block 
or a convolutional quantum code. Omitting the final Fourier 
transformation gates in the quantum circuit, we obtain an X- 

only stabilizer S' o = (I0\0). 

Expanding the matrix G\(D) as semi-infinite matrix, we 
get the following semi-infinite version of the stabilizer matrix 
G(D) of eq. ©: 



/ 911S2 
521 S2 



3l2<S'2 
3225*2 



ffl.ni 52 
32,711 52 
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3n52 31252 

32l52 32252 



31, ni 52 

32, ni 52 



3fci,i52 Sfci,252 • ■ • Sii,m& 
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This matrix indicates that we have to apply the inverse encod- 
ing circuit of the code C to every block of qudits corresponding 
to the submatrices gijS 2 . This first step corresponds to the 
leftmost boxes marked BC in the example of Fig. |4] The 
stabilizer matrix is now of the form 

G'(D) = Gi(D) <Z) P (1 0|0) = (Gi(D) $ 7|0). (9) 

This X-only generator matrix corresponds to a CSS code (see 
eq. ©) where 



H 2 (D) 



( Gi(D) 
\ 



\ 



Gi(D) J 

Using the algorithm of Sect. [Til] we obtain an inverse encod- 
ing circuit for the convolutional CSS code corresponding to 



3 



G\(D). This circuit has to be repeated r times if the identity 
matrix in eq. (O has rank r. The j-th copy of this quantum 
circuits acts on qudits j, j + r, j + 2r, . . . (see the blocks 
marked CCj in the example of Fig. |4j. Overall, we obtain 
the following result. 

Theorem 5: Let C be a quantum convolutional code which 
has been constructed using the product code construction 
described in Theorem |4] Denote the frame size with n and 
the constraint length with v. Then C has an encoding circuit 
whose depth is finite, i. e., does not depend on the length of 
the input stream. Furthermore, the depth of this circuit is upper 
bounded by poly(n, v). 

C. QCCs from products of cyclic codes 

In [8, Theorem 8], we have shown that product codes 
based on Reed-Solomon codes achieve the upper bound on 
the minimum distance of the resulting quantum code. Here 
we consider the following variant: 

Theorem 6: Let C be a cyclic code over ¥ q of composite 
length n = niri2 with n2\(q— 1). Furthermore, we assume that 
C can be decomposed as C = C\®C2 where C2 has generator 
polynomial g2(X) = n»=i (-^ — a *) wnere d—l< nz/2 and 
a is an 712-th root of unity in ¥ q . Then the code C is self- 
orthogonal and has dual distance d 1 - > min(<5i, d) where Si 
is the BCH bound of C^. 

Proof: The code C2 is generated by 

v»2-l 



Go 
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T,2(n 2 - 



\ 



,4-1 



v 2(d-l) 



,,(d-l)(n 2 -l) 



The inner product of row i and row j of G2 is 

ri2— 1 

as i + j ^ mod ri2. Hence C2 is self-orthogonal and so 
is C. The bound on the minimum distance follows from the 
two-dimensional BCH bound [1, p. 320]. ■ 
Starting with a generator matrix G = G\ ® G2 of a 
(permuted) cyclic code as in Theorem |6l we can construct 
convolutional quantum codes of CSS type. The semi-infinite 
generator matrix of the corresponding self-orthogonal convo- 
lutional code is formed by the copies of the generator matrix 
G which overlap in /J7i2 positions. For // = 2 we get 



/ (711G2 (/12G2 
321 G2 giiGi 



S11G2 
321 G2 



gl, II! G2 
52,rii G2 



gk 1 ,n 1 G2 

guG2 

<?22G2 



gi,n 1 G2 

52, ni G2 



gk 1 ,\G2 gfc 1 ,2G2 ■ ■ ■ gk 1 ,n 1 G2 



(10) 



The inner product of any two rows of this matrix is zero, 
as already G2 ■ G\ — 0. The dual distance of the dual of 
the convolution code defined by eq. dTOb is lower bounded 
by the dual distance of C, as any sequence in the dual of 
the convolutional code fulfills the parity checks given by the 
matrix G. Note that the encoder corresponding to the matrix 
([Tol l might be catastrophic. Then, in some cases, the minimal 
non-catastrophic encoder can have constraint length zero, i. <?., 
corresponds to a block code. 

VI. Example 

We illustrate the product construction and the corresponding 
encoding circuit using the five qubit code [5, 1,3] 2 and a 
classical convolutional code of rate R = 2/3. 
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Fig. 2. Inverse Encoding Circuit for the five qubit code. 

Using the inverse encoding circuit shown in Fig. |2] the 
stabilizer 
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of the five qubit code is transformed into 



S' = 



Note that the two Hadamard transforms on the fourth qubit 
cancel, but when omitting the final four Hadamard transfor- 
mations, we obtain an X-only stabilizer. 

In [10, Table 8.14] we find a nonsystematic rate R = 2/3 
convolutional code with memory v = 2 and free distance 
di iee = 3. An encoding matrix for that code is 

_ ( D + D 2 1 1 + D 2 



D + D 2 1 + D + D 2 



A minimal polynomial generator matrix for the dual code is 
given by 



1 + D + D 4 
H(D) = ( 1 + D 2 + D 3 + D 4 

l + D 2 + D A 



(12) 



If we apply our algorithm to the stabilizer code with X-only 
stabilizer matrix (X(D)\Z(D)) := (H(D)\0) we obtain the 
circuit shown in Fig. [3] The resulting transformed stabilizer is 
of the simple form (XII) 

Using the product construction, we take the tensor product 
of the stabilizer matrix S in eq. (|TTJ and the generator matrix 
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Fig. 3. Quantum circuit transforming the stabilizer (X(D)\Z(D)) :- 
(H(D)\0) into the simple form (XII). 



H(D) of the binary convolutional code in eq. ( fT2b . The 
stabilizer matrix has the form 



product = (H(D) ® S x | H(D) ® S z 



(13) 



where Sx and Sz denote the corresponding parts of S. 

Note that the circuit shown in Fig. [2] corresponds to a 
binary symplectic matrix T = T t T 2 , i.e., ST = S\ where T 2 
corresponds to the last four Hadamard gates. Replicating the 
circuit without these Hadamard gates three times as indicated 
in Fig. |U we get the matrix J3 £§) T\, where I3 denotes a 3 x 3 
identity matrix. Now the Z-part of the stabilizer is zero, and 
the X-part has the form 



1 + D + D A 
V D 2 + D 3 + D 4 
1 + D 2 +D i 



/I 
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1 0/ 



So replicating the circuit of Fig. [3] four times (but spread out 
to any every fifth qubit), we get an X-only stabilizer. The 
final four Hadamard gates in Fig. [4] transform it into a Z-only 
stabilizer. 

The structure of the whole encoding circuit is illustrated 
in Fig. |4] Only the first block is shown, but every quantum 
gate in the circuit has to be applied repeatedly, shifted by 
the corresponding number of qubits. Each block encodes 11 
qubits into 15. The inputs marked with bc^ correspond to the 
input of the i-th copy of the block code [5,1,3], the inputs 
of the four copies of the convolutional code are marked with 
cc^\ The boxes marked with BC correspond to the encoder 
for the block code in Fig. [2] the blocks CCj correspond to 
the encoding circuit for the convolutional code in Fig. [3] 

VII. Conclusions 

The problem of constructing quantum convolutional codes 
and their encoders was addressed. Using a CSS-type construc- 
tion, we derived new examples of QCCs of rate 2/4. For 
constraint lengths up to v = 6 we performed an exhaustive 
search of the search space, and for constraint lengths up to 11 
we employed a randomized search which found several good 
codes. Using a product code construction which takes as inputs 
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Fig. 4. Schematic inverse encoding circuit for the quantum convolutional 
code of rate _R = 11/15 obtained by the product code from the quantum 
block code [5, 1, 3]2 and a classical convolutional code with rate R = 2/3. 



a classical convolutional code on the one hand and a quantum 
block code on the other, it is possible to derive many examples 
of QCCs. We show that these codes all have the property that 
their encoder is of polynomial depth. We conjecture that any 
stabilizer QCC has a polynomial depth encoder. It seems that 
a more detailed study of the algorithm given in [9], which 
is based on iterative Smith normal form computation on the 
stabilizer matrix, would be required to resolve this question. 
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